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A two-locus genetic interaction between LPHN3 and 
11q predicts ADHD severity and long-term outcome 

MT Acosta 1 ' 4 , Jl Velez 1,4 , ML Bustamante 1 ' 2 ' 3 ' 4 , JZ Balog 1 , M Arcos-Burgos 1 and M Muenke 1 

The severity of attention-deficit/hyperactivity disorder (ADHD) symptoms is a major predictor of long-term ADHD outcome. 
To investigate if two-locus interactions might predict ADHD severity, we studied a sample of 1341 individuals from families 
clustering ADHD, using the Vanderbilt Assessment Scale for Parents. Latent class cluster analysis was used to construct 
symptom profiles and classify ADHD severity. Single nucleotide polymorphisms (SNPs) spanning ADHD-linked chromosomal 
regions on chromosomes 4, 5, 10, 11, 12 and 17 were genotyped. SNPs associated with ADHD severity were identified and 
potential two-locus genetic interactions were tested. We found that SNPs within the LPHN3gene interact with SNPs spanning the 
1 1 q region that contains DRD2 and NCAM1 not only to increase the risk of developing ADHD but also to increase ADHD severity. 
All these genes are identified to have a major role in shaping both brain development and function. These findings demonstrate 
that genetic interactions may predict the severity of ADHD, which in turn may predict long-term ADHD outcome. 
Translations! Psychiatry (2oW) 1, e17; doi:1 0.1 038/tp.201 1.14; published online 5 July 2011 



Introduction 

Attention-deficit and hyperactivity disorder (ADHD) is one of 
the most frequent psychiatric disorders in childhood, affecting 
10% of children and adolescents in the United States. 1 ADHD 
is characterized by elevated levels of inattention and/or 
hyperactive or impulsive behaviors that cause significant 
impairment 2 on a child's academic and social functioning. 3 

Several studies show that the severity of ADHD symptoms 
predicts long-term prognosis. 4-7 For example, Barry et al. 
reported that the strongest predictor of discrepancy between 
predicted achievement evaluated by measurements of 
intellectual function and actual academic performance is the 
severity of ADHD symptoms, defined as the index composed 
of averaged behavioral data obtained from the hyperactivity/ 
impulsivity checklist of DSM-IV, attention problems and 
hyperactivity scores from the Behavior Assessment System 
for Children-Parent Rating Scale, and from the total score on 
the Home Situations Questionnaire. 4 Similar results have 
been reported showing that the severity of inattention and 
hyperactivity negatively correlates with academic perfor- 
mance and has deleterious effects on other aspects related 
to quality of life, such as risk avoidance and satisfaction. 5 
Another study showed that ADHD severity was positively 
associated with the use of tobacco, alcohol and marijuana. 6 
In individuals without frank ADHD, higher levels of inattention 
have been found to have a negative effect on both academic 
performance and college adjustment. 7 

Due to this evidence, using a clinical approach that takes 
severity into account is relevant. Furthermore, the traditional 
approach using the diagnostic criteria for ADHD as defined in 



the DSM-IV discards those individuals who, despite the severity 
of some of the specific symptoms, do not fit the overall numeric 
criteria for the diagnosis of ADHD. However, the severity of their 
symptoms could have a marked functional impact and these 
individuals may benefit from interventions. 8-10 

In the etiology of ADHD, genetic factors are strongly 
implicated. 11 Recently, we identified a common variant 
of the Latrophilin 3 gene (LPHN3) that increases the risk of 
developing ADHD by 1.2-fold. The ADHD population attribu- 
table risk (PAR%) attributed to this gene is ~9%, suggesting 
that incidence of ADHD in the general population would be 
reduced by about 9% if it were possible to control the effect of 
the LPHN3 variant that confers susceptibility to ADHD. 12 
Functional studies revealed that LPHN3 is expressed in key 
brain regions related to attention and activity, and that its 
variants affect metabolism in neural circuits implicated in 
ADHD and are associated with response to stimulant 
medication. These findings were replicated in samples 
obtained from unrelated populations. 12 ' 13 In addition to 
increasing knowledge of the genetic basis of ADHD, the 
discovery of the molecular substrates of ADHD opens up new 
avenues for the exploration of targeted therapies. 

In related studies, we analyzed two-locus interactions 
models between LPHN3 variants on chromosome 4q and a 
susceptibility haplotype on 11q encompassing the NCAM1, 
TTC12, ANKK1 and DRD2 genes. The simultaneous pre- 
sence of both genetic risk variants increases the risk of ADHD 
by ~2.5 compared with the risk when none of these variants 
are present. 14 In addition to previous evidence linking these 
two regions with human behavior, 15-20 we found that this 
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interaction also explains differences in brain metabolism as 
observed through proton magnetic resonance spectroscopy 
( 1 H-MRS) data and pharmacogenetic response to stimulant 
medication better than the presence of either variant alone. 

In the present study, we aimed to evaluate the hypothesis 
that the interaction between ADHD susceptibility factors on 4q 
and 1 1q, as well as other genomic regions, not only predicts 
ADHD susceptibility but also the severity of the disease, as 
ADHD severity is the best indicator of long-term outcome. 
In this way, we additionally expand upon our previous findings 
to place them in a context where they may have a greater 
impact in terms of defining interventions for affected and 
subsyndromal individuals. 

Subjects and methods 

Subjects. We analyzed a total of 349 nuclear families 
consisting of a total of 1371 individuals. Participants 
were from a United States population, 4-65 years of age, 
ascertained from ADHD probands. Diagnosis of ADHD in 
children was established using the DSM-IV criteria; 2 in adults, 
the Conners Adult ADHD Rating Scale (CAARS) 21 was used. 
All participants were evaluated using the Vanderbilt 
Assessment Scale for Parents (VAS-P), 22 which includes all 
18 DSM-IV criteria for ADHD (questions 1-9 for inattention; 
questions 10-18 for hyperactivity/impulsivity), oppositional 
defiant disorder (ODD; questions 19-26), conduct disorder 
(CD; questions 27-40) and seven items from the Pediatric 
Behavior Scale 23 screening for anxiety and depression (A/D; 
questions 4^-47 in the VAS-P). Each question of the VAS-P 
is scored based on the frequency of presentation of the 
behavior on a 1 to 4 scale (1 = never, 2 = occasionally, 
3 = often and 4 = very often). 22 Complete information was 
available for 1341 individuals. More detailed information 
regarding clinical assessment and specific demographics of 
this sample is described elsewhere. 8 

Definition of severity. Latent class cluster analysis (LCCA) 
models containing 1-10 classes were fitted to the data using 
Latent GOLD 4.5 (Statistical Innovations, Belmont, MA, 
USA). Latent GOLD uses both EM and Newton-Raphson 
algorithms to find the maximum likelihood for each model 
after estimating model parameters. 24 In the LCC models, 
separate analyses were performed for each of the five VAS-P 
domains (inattention, hyperactivity, ODD, CD and A/D). The 
number of clusters was selected using a likelihood ratio test 
evaluating whether increases in likelihood (L 2 ) associated 
with increased latent classes justified their inclusion. 
Certainty of these clusters was assessed calculating 
P-values associated with L 2 values after running 500 
parametric bootstrap replicates. As covariates for all 
models, we used gender, ADHD diagnosis and age as 
categorical variables (children: 4-11 years; adolescents: 
12-17 years; adults >17 years). 12 

As implemented in Latent GOLD, individuals are assigned 
posterior membership probabilities for belonging to each 
cluster based on their symptom profiles. Individuals are 
additionally assigned to the cluster for which the posterior 
probability is highest. Based on this assignment, VAS-P 



profiles are obtained for each cluster. Using these profiles, 
a severity scale was derived as follows: consider the LCC 
results for a particular VAS-P domain, for example, inatten- 
tion, and let K be the number of clusters found, Q be the 
number of questions in that domain, n, be the number of 
individuals in that cluster (/= 1 ,2, . . . ,K) and N' P be the number 
of questions in cluster / with an average profile P. If N' P > Q/2, 
all n, individuals in cluster / will be classified as being severe for 
this domain, and not severe otherwise. In other words, 
individuals with severe symptoms will be those for whom the 
cluster they belong to, has at least Q/2 + 1 questions with a 
profile above 'Ocassionally'. From now on, individuals with 
severe symptoms will be treated as 'cases' and non-severely 
affected individuals as 'controls'. For the analyses, this 
derivation allows more precise descriptions of each indivi- 
dual's characteristics (profiles), descriptions that would 
be difficult to define when relying solely on the categorical 
DSM-IV diagnosis. 

Genotyping and genetic statistical analysis. DNA was 

prepared from peripheral-blood specimens and the lllumina 
genotyping service was used for genotyping, as previously 
described. 12 Out of the 1341 individuals who completed 
the VAS-P questionnaire, blood samples were available for 
1 181; genetic data was available for 813 individuals. A total 
of 369 markers homogenously distributed on genomic 
regions of chromosomes 4 (156 markers), 5 (48 markers), 
10 (7 markers), 11 (129 markers), 12 (15 markers) and 17 
(14 markers), as described previously, 12 were selected for 
analysis because of the presence of loci linked to ADHD in a 
genetic isolate from Colombia. 25 

Family-based association tests (FBATs) using a dominant 
model under the hypotheses of no linkage and no association 
were performed as implemented in Golden Helix's PBAT 
module (Golden Helix, Bozeman, MT, USA. Golden Helix 
PBAT Software http://www.goldenhelix.com). These hypoth- 
eses were selected because the linkage studies had been 
performed on a different population, and hence there was no 
previous knowledge about linkage or association in the United 
States population described here. As phenotypes, the VAS-P 
score, defined as the summation of the individual's responses 
in a particular domain, and the VAS-P severity scale, derived 
using LCCA, were used independently. The significance of 
markers was determined by using both the raw P-value and a 
permutation-based approach with B= 10000 permutations. 

Interactions effects on the severity of symptoms. We 

used the generalized Cochran-Mantel-Haenszel test as 
implemented in R, 26 to determine for all possible pairs of 
markers located on different chromosomes and significantly 
associated either with the VAS-P score or the severity of 
symptoms, whether the pair of single nucleotide 
polymorphisms (SNPs) are independent, conditionally on 
the levels of the severity phenotype (severe and not severe). 
This approach was stratified by domain, and using the 
definition of severity of symptoms previously described. A 
total of p pairs of markers were found to be of potential 
interest. By domain, the false discovery rate (FDR) 27 was 
used to correct by multiple testing. 
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To perform the Cochran-Mantel-Haenszel test for a 
particular domain we used the severity scale as the stratum 
indicator, such that a series of p pairs of contingency tables, 
each of dimension rx c, can be obtained. In our context, p 
(p = 587) is the number of possible pairs of markers located 
within different chromosomes, r (r=3) is the number of 
possible genotypes in marker 1 and c (c=3) the number of 



possible genotypes in marker 2. To avoid any bias inherent to 
family data, we selected only those families with (i) discordant 
parents and concordant children affected, (ii) two discordant 
siblings and unaffected parents, (iii) unaffected parents and 
several discordant children and (iv) discordant parents and 
discordant siblings. Families falling in the third case were 
further examined so only one affected and one unaffected 



Table 1 Summary of clusters derived by latent class analyses 



Questions Criteria? 



Cluster Characteristic and affection status for 
broader phenotype 



Predominant gender 
and age classes 



Severity of 
symptoms 



n (%) 



Inattention 
1-9 >5 



Hyperactivity/impulsivity 
10-18 >5 



ODD 
19-26 



>3 



Anxiety and depression 
41-47 >2 



Minimal symptoms, unaffected for ADHD 
Few symptoms, unaffected for ADHD 
Mostly inattentive, high presence of symptoms, 
affected for ADHD 

Mostly inattentive, higher presence of 
symptoms, affected for ADHD 
Fewer symptoms, affected for ADHD 
Fewer symptoms, affected for ADHD 



Minimal H/l symptoms, most individuals 

unaffected for ADHD 

Minimal H/l symptoms, most individuals 

unaffected for ADHD 

Few symptoms, affected for ADHD 

High presence of symptoms except for Q12, 

most individuals affected for ADHD 

Few symptoms for Q1 0-1 4 and fewer for 

Q15-Q18, individuals are mostly affected 

Few symptoms for Q10 and Q1 1 , but minimal 

symptoms for the rest; ADHD affection status 

is equally present 

Mostly H/l individuals, all affected for ADHD; 

higher presence of symptoms 

Mostly H/l individuals affected for ADHD; 

high presence of symptoms 

Few-to-minimal symptoms except in Q1 5; ADHD 

affection status is equally present 



Minimal symptoms, most individuals unaffected 
for ADHD 

Few symptoms, most individuals affected for 
ADHD 

High presence of symptoms in Q19-25 and few 

for Q26; most individuals affected for ADHD 

Minimal symptoms in all questions; individuals 

are mostly unaffected for ADHD 

Minimal symptoms in all questions; ADHD 

affection status is equally present 

Few symptoms in Q20 and Q24; ADHD affection 

status is equally present 

Higher presence of symptoms in all questions; 

most individuals affected for ADHD 



Minimal symptoms in all questions; individuals 

are mostly unaffected for ADHD 

Fewer symptoms; mostly affected individuals 

Few symptoms; most individuals unaffected for 

ADHD 

Fewer symptoms in Q41 , Q42 and Q47 with few 
on the rest; most individuals affected for ADHD 
High presence of symptoms; most individuals 
affected for ADHD 

Higher presence of symptoms in all questions; 
most individuals affected for ADHD 



Females-adults 
Males-adults 
Males-children and 
adolescents (few adults) 
Males-children and 
adolescents (few adults) 
Males-adults 
Males-children 


Not severe 
Not severe 
Severe 

Severe 

Severe 
Severe 


370 (27.6) 
313(23.4) 
271 (20.2) 

182 (13.6) 

151 (11.3) 
53 (4.0) 


Females-adults 


Not severe 


418(31.2) 


Males-adults 


Not severe 


307 (22.9) 


Males-adolescents 
Males-children and 
adolescents (few Adults) 
Females-adults 


Severe 
Severe 

Severe 


165 (12.3) 
165(12.3) 

101 (7.5) 


Malp^-pll aripQ 

IVICMCO GUI uUCO 


Nnt qpvptp 


58 (4.3) 


Males-children 


Severe 


53 (4.0) 


Males-children 


Severe 


48 (3.6) 


Males-children 


Not severe 


25(1.9) 


Female-adults 


Not severe 


365 (27.2) 


Male-all ages 


Not severe 


293 (21.9) 


Males-children 


Severe 


182(13.6) 


Males-adults 


Not severe 


161 (12) 


Males-children and 
adolescents (few adults) 
Females-adults 


Not severe 
Not severe 


149 (11.1) 
143(10.7) 


Males-children and 
adolescents (few adults) 


Severe 


47 (3.5) 


Males-adults 


Not severe 


375 (28.0) 


Females-adults 
Males-children 


Not severe 
Not severe 


376 (28.1) 
255 (19) 


Males-children 


Not severe 


157(11.7) 


Females-adults 


Severe 


132 (9.9) 


Females-adults 


Severe 


45 (3.4) 



Abbreviations: ADHD, attention deficit hyperactivity disorder; adolescents, 1 2-1 7 years; adults, > 1 7 years; children, 4-1 1 years; H/l, hyperactivity/impulsivity; ODD, 
oppositional defiant disorder; Q, question. 

a Number of questions being marked as 'often' or 'very often' in the Vanderbilt Assessment Scale for Parents (VAS-P) questionnaire. 
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child was selected. This selection was performed matching by 
gender and age range. Parents were excluded from families 
falling in the fourth case. At the end of the process, only one 
discordant pair was selected within those families that fulfilled 
the inclusion criteria. 



all of the independent variation conferred by the spectrum of 
symptoms that define ADHD and its comorbidities. Further, 
using this categorization for each domain, we performed 
association and interaction analyses for the severity 
phenotype. 



Results 

Severity of ADHD symptoms. In our maximally expanded 
sample (n=1341), we identified six significant independent 
and mutually exclusive clusters for the VAS-P inattention 
domain, nine for hyperactivity/impulsivity, seven for ODD and 
six for the A/D (Table 1 and Figure 1), all of which were 
consistent with findings from previous studies. 8,12,28-31 
Because of small variability in the predictors in the CD 
domain, LCCA models could not be fitted to this data. 
Individuals were considered as having severe symptoms 
when belonging to clusters 3, 4, 5 and 6 in the inattention 
domain (n = 647, 49%); to clusters 3, 4, 5, 7 and 8 in the 
hyperactivity/impulsivity domain (n = 532, 39.7%); to clusters 
3 and 7 in the ODD domain (n = 229, 17%) and to clusters 5 
and 6 in the A/D domain (n=177, 13.2%) (Figure 1). Under 
this categorization, severe versus not severe status involves 



FBAT. Once the severity of symptoms for each of the four 
VAS-P domains was identified, an FBAT was performed, 
taking into account the complete family structure and the 
severity status (severe and not severe) as disease indicator 
(outcome). In complementary analyses, the overall VAS-P 
score for each dimension was also used as the outcome for 
the FBAT as a general indicator of severity. 

Table 2 presents the FBAT results for both the LCCA- 
derived severity status and the overall VAS-P score for each 
domain. After permutation, a total of 38 genetic markers were 
found to be significantly associated either with the VAS-P 
score or the severity of symptoms by dimension. 

In the inattention domain, a total of eight markers were 
found to be associated with either the LCCA-derived severity 
status or the VAS-P score. When considering the VAS-P 
score, significant associated markers were identified within 
LPHN3(4q), NDFIP1 (5q), and DRD2 (1 1q) and TTC12{\ 1q); 





— Cluster 1 — Cluster 4 Cluster 7 

— Cluster 2 Cluster 5 Cluster a 

— Cluster 3 — Cluster© — Clusters 



Figure 1 Profile plots derived using latent class cluster analysis applied to attention-deficit and hyperactivity disorder symptoms as measured by the Vanderbilt 
Assessment Scale for Parents (VAS-P) questionnaire. For each domain, VAS-P symptoms profiles within clusters are shown in a scale from 1 to 4 (1 = never, 
2 = occasionally, 3 = often and 4 = very often). Demographic characteristics for each cluster are shown in a scale from 0 to 1 , representing the proportion of individuals from 
the population with such characteristics, (a) Inattention (Q1-Q9): Q1 : careless, inattentive; Q2: sustains attention poorly; Q3: appears to not listen; Q4: poor follow through; Q5: 
disorganized; Q6: avoids/dislikes sustained mental effort; Q7: loses needed objects; Q8: easily distracted; and Q9: often forgetful, (b) Hyperactivity/impulsivity items ( Q10- 
Q18): Q10: fidgets or squirms; Q11 : cannot stay seated; Q12 restless; Q13: loud; noisy; Q14: always 'on the go'; Q15: talks excessively; Q16: blurts out; Q17: impatient; and 
Q18: intrusive, (c) Oppositional defiant disorder (ODD; Q19-Q26): Q19: loses temper; Q20: argues with adults; Q21: defies adults' rules; Q22: annoys others; Q23: shifts 
blames to others; Q24: touchy; Q25: angry/resentful; and Q26: vindictive, (d) Anxiety and depression (Q41-Q47): Q41: fearful, worried; Q42: fear of making mistakes; Q43: 
feels useless; Q44: blames self; Q45; feels unloved; Q46: sad; and Q47: embarrassed. 
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Table 2 Raw and permuted P-values from family-based association tests using the VAS-P score and the severity of symptoms derived using latent class analyses, as 
qualitative and quantitative traits, respectively 



Marker 


Gene 


Allele 


Frequency 


Effei 


Chromosome 4 










rs1 947275 


LPHN3 


C 


0.814 


+ 


rs21 32074 


LPHN3 


G 


0.618 


+ 


Ib'fOO^OUU 


/_PHA/3 


ri 


U.DUO 




rs1 31 24636 


LPHN3 


G 


0.046 




rs4860091 


LPHN3 


T 


0.386 


+ 


rs335322 


LPHN3 


G 


0.575 


+ 


rs1 001 5239 


LPHN3 


A 


0.577 


+ 


rs351 06420 


LPHN3 


G 


0.985 


+ 


rs1 2646895 


LPHN3 


G 


0.508 


+ 


ib i o i uy^u 


/_PHA/3 




U.UD I 




rs1 86750 


LPHN3 


A 


0.243 


+ 


rs734644 


LPHN3 


T 


0.271 


+ 


rs6813183 


LPHN3 


G 


0.300 


+ 


rs6551670 


LPHN3 


A 


0.300 


+ 


rs6551669 


LPHN3 


C 


0.300 


+ 


Chromosome 5 










rs249637 


NDFIP1 


G 


0.074 


- 


rs37022 


SLC6A3 


T 


0.831 


+ 


rs250682 


SLC6A3 


C 


0.798 


+ 


rs6596271 


SLC25A48 


G 


0.042 


- 


rs1 64080 


GNPDA1 


C 


0.526 


+ 


Chromosome 11 








rs1 0891 551 


DRD2 


A 


0.122 


+ 


rs719804 


TTC12 


G 


0.223 


- 


rs4938006 


NCAM1/ 


G 


0.109 


+ 




TTC12 








rs1 2799083 


DRD2 




U.UO I 




rs4245148 


DRD2 


T 


0.122 


+ 


rs1 759601 7 


NCAM1 


T 


0.024 


- 


rs1381246 


BSX 


p 


U.'fO 1 




rs652285 


NCAM1 


T 


0.057 


+ 


rs675646 


NCAM1 


C 


0.071 


+ 


rs649568 


NCAM1 


T 


0.044 


+ 


rs702966 


PHRF1 


C 


0.276 




rs2574829 


NCAM1 


G 


0.394 


+ 


rs1 1214521 


NCAM1 


G 


0.035 




rs1 1214505 


NCAM1 


G 


0.082 




rs1 055076 


TTC12 


A 


0.175 


+ 


rs1 2222469 


NCAM1 


A 


0.035 


+ 


Chromosome 12 








rs1 77221 34 


TPH2 


G 


0.029 




Chromosome 17 








rs7208257 


ARRB2 


C 


0.055 


+ 



Score 



Severity 



Inattention 



0.0081 9 a 



H/l 



ODD 



0.001 7 a 
0.0031 b 
0.0414 b 



CD 



0.0054 b 

0.007 b 

0.0098 b 



A/D Inattention 



H/l 



ODD 



A/D 



0.01 54 b 
0.001 9 b 
0.0368 b 
0.021 6 b 



0.01 17 b 
0.01 14 b 
0.01 05 b 
0.01 05 b 



0.021 6 b 



0.0253 b 
0.01 17 b 



0.0052 a 
0.0082 b 



0.0052 b 



0.034 b 0.01 06 a 

0.01 07 b 
0.0251 b 



0.0191° 
0.0057 b 



0.0097 b 
0.01 07 b 
0.0171 b 
0.0121 b 



0.0082 b 
0.0807 b 



0.0084 b 



0.0226 
0.0292 
0.0343 



0.0015° 



0.0359 a 



0.0214 b 
0.01 05 a 
0.01 62 b 



0.01 83 b 



0.0085 b 



0.01 64 b 



Abbreviations: A/D, anxiety and depression; CD, conduct disorder; H/l, hyperactivity/impulsivity; ODD, oppositional defiant disorder; VAS-P, Vanderbilt Assessment 
Scale for Parents. Nonsignificant P-values after B = 1 0,000 permutations are shown in italic. No correction for the number of questions interrogated was performed in 
genetic association analyses. 
a Permutated P-value <0.01. 
b Permutated P-value <0.05. 

c Effect on the severity towards opposite direction to the one described in the 'effect' column. 



for the severity scale, all significantly associated markers are 
within NCAM1 (11q). 

For hyperactivity/impulsivity, two markers within SLC6A3 
(5p) were found to be associated with the VAS-P score and 
three markers within NCAM1 were associated with the 
severity of symptoms. 

For ODD, four markers were found to be associated with the 
VAS-P score and four with the severity of symptoms. For the 
former, three of the four markers are within LPHN3an6 one is 
within NCAM1. For the severity of symptoms, all significant 
markers are within LPHN3. 

For CD, six markers were found to be associated with the 
VAS-P score, three within LPHN3 and three within NCAM1. 

In the A/D domain, a total of eight markers were found to be 
associated with the VAS-P score and six with the severity of 



symptoms. Among the markers found to be associated with 
the VAS-P score, four are within LPHN3, two are within 
NCAM1, one is within GNPDA 1 and one is within TPH2 (12q). 
Regarding severity of symptoms, a total of six markers were 
found to be associated, four within NCAM1, one within 
GNPDA1, and one at an intergenic region between NCAM1 
and TTC1 2 on chromosome 11. 

LPHN3x 11q interaction. By domain, a total of 567 unique 
pairs of markers located at different chromosomes were 
tested for SNP x SNP interactions. After breaking up the 
family data using the strategy previously mentioned, the 
number of severe (cases)/not severe (controls) pairs of 
individuals available for analysis was 184 for the inattention 
domain, 169 for hyperactivity/impulsivity, 104 for ODD and 
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a 

Ice 



u 



AA 



AC CC AA AC 

rs1 759601 7 (NCAM1) 



CC 



£2 

|CC 



CT 



W 



AA 



AC CC AA AC 

rs1 2799083 (DRD2) 



CC 



i= GG 



I AG 



n AA 



U 



AA 



AC 



CC AA 
rs620291 (NCAM1) 



AC 



CC 



CT 



W 



AA 



AG GG AA AG 

rs1 121 4505 (NCAM1) 



GG 



CT 



CO 



§CC 

r, AC 



£ AA 



U 



0.40 
0.35 
0.30 
0.25 
0.20 
0.15 
0.10 
0.05 
0.00 



AA 



AC CC AA AC 

rs620291 (NCAM1) 



CC 




rrO.8 

0.7 
0.6 
0.5 
0.4 
0.3 
0.2 

|o, 

Lr°-° 



AA 



AC CC AA AC 

rs4938006 (NCAM1) 



CC 



76 for A/D. Figure 2 presents the results of the Cochran- 
Mantel-Haenszel tests after correcting by multiple testing 
using the FDR. Regardless of the VAS-P domain being 
considered, statistically significant SNP x SNP interactions 
include markers within LPHN3 and within the region 



Figure 2 Genotypic frequency distribution for pairs of single nucleotide 
polymorphisms (SNPs) contributing to interaction effects and assessed by the 
Cochran-Mantel-Haenszel test. Significant interaction effects involve markers in 
the LPHN3 gene and a region in 1 1 q that harbors the NCAM1 and DRD2 genes. 
Epistatic effects are depicted by changes in color that represent significant 
differences in genotypic distribution among severe (cases) vs not severe (controls) 
individuals. In there, the genotypes for one marker are held fixed whereas 
genotypes on the other marker vary. Inattention: (a) markers rs1 947275 harbored in 
LPHN3 and rs1 759601 7, in NCAM1, contribute to the severity of symptoms in this 
domain (11^=33.163, FDR-corrected P-value <0.001); (b) markers rs1947275 
harbored in LPHN3 and rs1 2799083, in DRD2, produce a significant interacting 
effect contributing to the severity of symptoms (h/P = 28.456, FDR-corrected 
P-value < 0.005). H/l: (c) markers rs351 06420, in LPHN3, and rs620291, in 
NCAM1, produce an interacting effect contributing to the severity of symptoms 
(^ = 20.497, FDR-corrected* P-value <0.05). ODD: (d) markers rs995447, in 
LPHN3, and rs1 1214505, in NCAM1, interact to modify the severity of symptoms 
(^ = 41 .379, FDR-corrected P-value < 0.0001); (e) rs734644, in LPHN3, and 
rs4938006 produce an epistatic effect contributing to the severity of symptoms 
(^ = 26.795, FDR-corrected P-value <0.01). A/D: (f) markers rs1510920, in 
LPHN3, and rs4938006 localized closest to NCAM1 in an intragenic region of 
chromosome 1 1q, interact to modify the severity of symptoms (Vtf 2 = 41 .379, FDR- 
corrected P-value < 0.0001). *FDR-corrected P-values calculated on the basis 
of 567 independent tests, corresponding to the maximum number of SNP pairs 
for each of the four domains of the VAS-P questionnaire from which the severity of 
symptoms was derived. Abbreviations as in Figure 1. 

containing NCAM1, DRD2 and TTC12 genes (Table 2 and 
Figure 2). 

For the inattention domain, the interaction between 
rs1 947275 and rs1 759601 7 was found to be statistically 
significant (IVp = 33.163, P F DR-corrected = 0.00067). In this 
case, comparison of the genotypes distribution for these 
markers between cases and controls showed that decreasing 
the number of copies of the T allele in LPHN3 (from CT to CC 
in rs1 947275) yields a reduction of the number double 
homozygous individuals for the C allele (marker 
rs1 759601 7) (Figure 2). 

For the hyperactivity/impulsivity domain, in which the 
interaction between rs351 06420 and rs620291 was signifi- 
cant (^ = 20.497, P F DR-corrected = 0.02766), similar results 
were obtained. For instance, when the G allele is present in 
LPHN3 and the allele C in NCAM1 (while keeping the other 
fixed), the number of cases with a particular genotype 
decreases when compared with the controls (Figure 2). 

For ODD, two interactions were found to be statistically 
significant, one defined by the markers rs995447 and 
rs1 121 4505 (^ = 41.379, P FD R-corrected = 0.00008), and 
one by the markers rs734644 and rs620291 (l\/P = 26.795, 
PFDR-corrected = 0.001 96). In the former, the presence of one 
copy of the G allele in NCAM1 yields a reduction in the number 
of cases with two copies of the C allele in LPHN3 and an 
increment of those with only one copy of it. In the latter, 
regardless the number of copies of the T allele in LPHN3, 
the number of cases with one or more copies of the C 
allele in NCAM1 is greater than the number of controls 
(Figure 2). 

In the A/D domain, only the interaction between rs1510920 
and rs4938006 was statistically significant (IVP = 23.973, 
PFDR-corrected = 0.00655). For instance, when the one copy 
of the C allele is present in NCAM1, it yields a reduction in the 
number of cases with one or more copies of the A allele in 
LPHN3 (Figure 2). 
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Discussion 

We found that several genetic markers interact to increase the 
severity of ADHD symptoms in our population. These 
markers, distributed along five chromosomal regions that 
have previously been linked to ADHD in a genetic isolate 
from Colombia, 25 are harbored in genes that have a potential 
role in the neurobiology of ADHD. In our sample, two 
genetic interactions, LPHN3x NCAM1 and LPHN3x DRD2, 
modified the severity of symptoms in ADHD. The former 
modifies the inattention, hyperactivity/impulsivity, ODD and 
A/D symptoms in ADHD, the latter acts on the inattention 
domain only. 

Severity was estimated in our sample by two methods, one 
using the overall VAS-P score, and the second quantifying the 
magnitude of inattention, hyperactivity/impulsivity, CD, ODD 
and A/D symptoms using LCC. By using this methodology, we 
could include in the analysis all the clinical information from 
the individuals regardless of their categorically diagnosed 
status. This approach may have a more direct translation to 
the clinical context, where severity of symptoms is a major 
predictor of functional outcome. 

Of special interest in our findings is the association and 
interactions found in regions that harbor genes previously 
described as being associated with ADHD, such as /.PHA/3 12 
and DRD2, 20 as well as the novel interaction between LPHN3 
and NCAM1 recently described by our group. 14 These genes 
are expressed in areas of the brain that have highly related 
roles, which makes them relevant candidates to explain the 
type and severity of ADHD symptoms. 

LPHN3 is a novel finding in association with ADHD. Our 
group has recently described the association between a 
specific LPHN3 haplotype and the increased risk of ADHD, 
particularly related to a higher risk of behavioral problems 
such as conduct disorder and substance use disorder. 12 In 
this study, we found that the interaction of three genetic areas 
highly related to neurodevelopmental aspects as well as 
function of the reward system predicts severity of ADHD 
symptoms. This is a very interesting observation that 
deserves further replication in a bigger sample. 

The region located in chromosome 11q23 includes the 
NCAM1, TTC12, ANKK1 and DRD2 genes, all of which are 
functionally linked to dopamine in the brain. Many association 
studies of DRD2 and substance dependence, including 
alcohol dependence and drug dependence have been 
reported. 32,33 

NCAM1 is expressed in all cell types in the nervous system, 
where it promotes neuron-neuron and neuron-glia adhesion. 
It is also involved in signal transduction and critical develop- 
mental events such as migration and proliferation, neurite 
outgrowth and fasciculation, synaptogenesis and synaptic 
plasticity 34,35 It has important roles in the development of the 
pre-frontal cortex 19 as well as in the regulation of dopaminer- 
gic pathways 15 and response to stress. 18 

Our previous results show that the association of the marker 
rs65511665 (harbored in the LPHN3 gene) is stronger in 
families from the Paisa population, a genetic isolate in 
Colombia that has a high prevalence of ADHD and a subtype 
highly co-morbid with disruptive behaviors. 12 ' 36 Furthermore, 
the same marker is associated with response to stimulant 



medication 12 ' 36 and, more recently, an interaction between 
this marker and an haplotype in 1 1 q increases the risk of 
developing ADHD ~ 2.5-fold. 14 The present study expands 
previous work by adding an additional clinically relevant 
distinction. The role of NCAM1 in ADHD has not been 
extensively evaluated. However, according to our current 
findings, this gene and the genetic markers in its vicinity 
deserve further evaluation in the context of the implications 
that it may have in the clinical expression of ADHD symptoms. 
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